THE DISTANCE FUNCTION FROM THE BOUNDARY 
IN A MINKOWSKI SPACE 



GRAZIANO CRASTA AND ANNALISA MALUSA 



Abstract. Let the space R" be endowed with a Minkowski structure M (that 
is AI : R" — > [0, +00) is the gauge function of a compact convex set having the 
origin as an interior point, and with boundary of class C^), and let d^^ (x, y) be 
the (asymmetric) distance associated to M. Given an open domain H C M" of 
class C^, let dfi{x) := inf{d^^ {x,y); y G dfl} be the Minkowski distance of a 
point X £ Q from the boundary of Q. We prove that a suitable extension of 
to M" (which plays the role of a signed Minkowski distance to 9f2) is of class 
in a tubular neighborhood of dfl, and that dfj is of class outside the cut 
locus of dfl (that is the closure of the set of points of non— differentiability of dQ 
in n). In addition, we prove that the cut locus of dil has Lebesgue measure 
zero, and that f2 can be decomposed, up to this set of vanishing measure, 
into geodesies starting from dQ and going into Q along the normal direction 
(with respect to the Minkowski distance). We compute explicitly the Jacobian 
determinant of the change of variables that associates to every point x £ Q 
outside the cut locus the pair {p{x),dQ{x)), where p{x) denotes the (unique) 
projection of x on dfl, and we apply these techniques to the analysis of PDEs 
of Monge-Kantorovich type arising from problems in optimal transportation 
theory and shape optimization. 



1. Introduction 

In recent years the study of the distance function from the boundary has at- 
tracted the attention of many researchers coming from different areas of mathemat- 
ical analysis. We mention, among other papers, in the framework of general 
first order Hamilton- Jacobi equations and Finsler geometry, iQ and in the 
case of Riemannian manifolds, and [pTf for the point of view of non-smooth 
analysis in Hilbert spaces, ||l3| for results applied to the theory of Sobolev spaces, 
and for applications to causality theory. 

The results of the present paper are mainly motivated by their applications to the 
analysis of PDEs. Let fl C K", n > 2, be a smooth open domain (i.e., a nonempty 
open bounded connected subset of M", with sufficiently smooth boundary), and let 
H : X R" ^ R be a smooth hamiltonian, such that for every x € the sublevel 
K{x) := {p G R"; H{x,p) < 1} is a compact convex set having the origin as an 
interior point, with boundary and with strictly positive principal curvatures at 
any point {K{x) G for short). It is well known that the function 

F{x,v) := ma.x{{v, p) ■,p e K{x)}, x eU, w G M" 
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is a Finsler structure on Q (see that is, a smooth function defined in f2 x (R" \ 
{0}), positively 1-homogeneous with respect to v, and such that the Hessian matrix 
dl^^_F'^ is positive definite at every point of x (R" \ {0}). 

For every pair of points a;,?/ G O, we can define the geodesic distance from y to 
X as ^ 

:=inf|^ F{m.m)dt- ^ e S^^y 

where Sx,y is the set of all Lipschitz curves ^: [0, 1] — > such that ^(0) — y and 
^(1) — X. It is well known that the distance from the boundary of f2, defined by 

(1.1) dn{x) -.^ inf L{x,y), xefl 

is the unique viscosity solution of the Hamilton- Jacobi equation 

(^2^ (Hix,Duix))^l inn, 

I M = on dfl 



(see |16|). Moreover, for every x E Q the infimum in (1.1) is achieved at a point 



y £ dfl that can be joined to x by a geodesic going into and starting from y along 



the "normal" direction to dfl (see [15 



It is worth to remark that, in the simple example H{x,p) = \p\, then F{x,v) 



\v\, (R",F) is the standard Euclidean space, and (1.2) is the eikonal equation 



Hence da is the Euclidean distance of x from dQ, that is 
dnix) — inf |a; — i/j , x £ fl . 

Another case, that will be the one considered in our analysis, concerns au- 
tonomous Hamiltonians of the form H{x,p) = p{p), where p{p) :— iid{t > 0; p £ 
tK} is the gauge function of a fixed compact convex set K £ C\. In this case the 
function F{x,v) coincides with the gauge function p°(f) of the polar set K'^ of K, 
and (R",F) is a Minkowski space (see ^, Chap. 14]). The function 

dnix) = inf p°{x ~ y) , x £ fl , 

which is the unique viscosity solution of 

{p{Du) = l inn, 
I u = on dVl , 

will be called the Minkowski distance from dfl. This non-symmetric distance func- 
tion from the boundary is exactly the object we shall deal with in this paper. 

Our results can be divided into two main groups, the first one devoted to the 
differentiability properties of da , and the other one on the regularity of the closure 
E of the set S of those points in where do is not differentiable. 

Concerning the differentiability of do,, we consider a signed distance d^, defined 
in the whole R", which extends do,, and we prove that d^ is of class in a tubular 



neighborhood of dn (see Theorem 4.16). This result provides a good definition of 



Dda and D^d^ on dn. For every y £ dn, Ddu{y) is proportional to the inward 



normal v{y) to dD. at y (see Lemma 4.3). Moreover, due to the Minkowskian 



structure of the space, for every x £ VL and for every projection a^o of x on dVL the 
geodesic jointing xq with a; is a segment starting from xq and going into SI along 



the "normal" direction, which is Dp{v{xo)) (see Proposition 4.4). Thanks to these 
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regularity results, for every y 6 we are able to define and study a linear mapping 
acting on the tangent space to dO, at y which plays the role of the Weingarten map 
in Riemannian geometry, that is it has real eigenvalues which are normal curvatures 
with respect to the Minkowski distance. 

For what concerns the singular set, a well known regularity result is the rectifi- 
ability of S. Namely, since dji is a locally semiconcave function on VL (see Q ^), 
then E is C^-rectifiable (see ) , that is it can be covered by a countable family of 
embedded manifolds of dimension rt — 1 , with the exception of a set of vanishing 
ffn-i ijieasure (here W denotes the s-dimensional Hausdorff measure). In general, 
this is the best regularity result that we can expect for E, and remains valid even 
without any smoothness assumption on dQ. 

On the other hand, the set E may behave badly. In it is exhibited a domain 

C M^, of class C^'^, such that E has positive Lebesgue measure. (In the same 
paper the authors prove a rectifiability result in a Riemannian setting.) Li and 
Nirenberg in |15 have shown that, if fl is of class C^'^ and F is of class C°° , then 
7i"~^(E) is finite, hence its Hausdorff dimension does not exceed n — 1. For related 
results in Finsler spaces see also [|[ §6] and 

Our result in this direction is the following (see Theorem 3.11): If il and 
F{x,p) — p^{p) are of class C^'", for some a e [0, 1], then the Hausdorff dimension 
of E does not exceed n — a. In the case a = (that is, Q. and of class C^), we 
show that E has vanishing Lebesgue measure, and that dn is of class on SI \ E 
(see Corollary |6 . 9| and Theorem 3.1C). Finally, the set S1\E can be decomposed into 
geodesies (segments) starting from and going into along the normal direction. 
We compute explicitly the Jacobian determinant of the change of variables which 
associates to every point a; £ \ E the pair (p(a;), dQ_{x))^ where p{x) denotes the 
(unique) projection of x on and we are able to perform a change of variables 
in multiple integrals in U, (see Theorem 7.1 ). 

As an application of these results, we consider the following system of PDEs of 
Monge-Kantorovich type: 



(1.3) 



-d\v{vDp{Du)) = f in f}, 
p{Du) < 1 in ri, 

p(Du) = 1 in{v> 0}, 



where the source / > is a continuous function in fJ, complemented with the 
conditions 



(1.4) 



f M > 0, w > in f], 
I M = on dVt. 



The first equation in (1.3) has to be understood in the sense of distributions, whereas 
w is a viscosity solution to the Hamilton-Jacobi equation p{Du) = 1 in the set 
{v > 0}. We look for a solution {u,v) to (L3)-(L4) in the class of continuous and 
non-negative functions. 

This system of PDEs arises in problems of shape optimization (see Q). In the 
case p{^) — 1^1, ( |l.3|)- (pr^) describes the stationary solutions of models in granular 
matter theory (see |d]). 

In section |^, using a change of variable formula, we shall explicitly construct 
a non- negative and continuous function Vf (defined in (7/7)) such that the pair 
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{dn,Vf) is a solution to (1. 

2. Notation and Preliminaries 

2.1. Basic notation. The standard scalar product of two vectors x,y ^ R" is 
denoted by {x, y), and \x\ denotes the Euclidean norm of x e M". By S*"^^ we 
denote the set of unit vectors of M", and by Aik the set oi kxk square matrices. We 
shall denote by (ei, . . . , e„) the standard basis of R". Given two points x,y G M", 
[x,y] will denote the closed segment joining x to y, while {x,y) will denote the 
same segment without the endpoints. li A, B C M", x G R" and t G R, we define 
A + x = {a + x; a€ A},tA^ {ta: a e A} and A + B =^ {a + b; a € A, b € B}. 

As is customary, Br{xo) and Br{xo) are respectively the open and the closed 
ball centered at xq and with radius r > 0. Given two vectors v,w £ M", the symbol 
V 1^ w will denote their tensor product, that is, the linear application from M" to 
M" defined by (ti (g) w){x) = v {w, x). 

We shall denote by C^{A) and H^{A), respectively, the Lebesgue measure and 
the s-dimensional Hausdorff measure of a set A C R" . 

Given A C R", we shall denote by Lip(A), C{A) and C'=(^), k e N, the sets 
of functions u : A — > R that are respectively Lipschitz continuous, continuous and 
fc-times continuously differentiable in A. Moreover, C°°{A) will denote the set of 
functions of class C''{A) for every k £ N, while C'^'"(A) will be the set of functions 
of class C'' (A) with Holder continuous k-th partial derivatives with exponent a G 
[0,1]. 

A bounded open set A C R" (or, equivalently, its closure A or its boundary dA) 
is of class C'^, fc G N, if for every point xq G dA there exists a ball B — Br{xo) and 
a one-to-one mapping -tp: B ^ D such that -tjj G C^{B), -ip'^ G C^{D), ^j{Br\A) C 
{x G R"; Xn > 0}, ip{B n OA) C {x G R"; a;„ = 0}. If the maps 'ip and ^p-^ are 
of class C°° or C'^''" (fc G N, a G [0, 1]), then A is said to be of class or C'^'" 
respectively. 

2.2. Differential geometry. We recall briefly some elementary facts from differ- 
ential geometry of hypersurfaces of class (see e.g. |^). Let A c R" be a 
bounded open set of class C^. For every x G OA, we denote respectively by j/(x) 
and TxA the unique inward unit normal vector and the tangent space of dA at x. 
The map i': dA S"'~^ is called the spherical image map (or Gauss map). It 
is of class and, for every x G dA, its differential di/a, maps the tangent space 
TxA into itself. The linear map := —diy^ ■ T^A — > T^A is called the Weingarten 
map. The bilinear form defined on T^A by Sx{v, w) — {Lx v, w), v,w € TxA, is the 
second fundamental form of dA at x. The geometric meaning of the Weingarten 
map is the following. For every v G TxA with unit norm, Sx{v,v) is equal to the 

normal curvature of dA at x in the direction v, that is, Sx{v,v) = (^£,{0), v{x)^, 

where ^(i) is any parameterized curve in dA such that (^(0) = x and ^(0) = v. 
The eigenvalues ki{x), . . . ,Kn~i{x) of the Weingarten map Lx are, by definition, 
the principal curvatures of dA at x. The corresponding eigenvectors are called the 
principal directions of dA at x. It is readily shown that every Ki{x) is the normal 
curvature of dA at a; in the direction of the corresponding eigenvector. From the 
regularity assumption on the manifold dA, it follows that the principal curvatures 
of dA are continuous functions on dA. 



)-(1.4), extending a result proved in M in the case 
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2.3. Convex geometry. By /Cq we denote the class of nonempty, compact, convex 
subsets of M" with the origin as an interior point. We shall briefly refer to the 
elements of /Cq as convex bodies. The polar body of a convex body K G /Cq is 
defined by 

= {pe M"; {p, x) <iyxeK}. 

We recah that, if K e /C^, then e /C^ and = [K^f = K (see [H 

Thm. 1.6.1]). 

Given K G /Cq we define its gauge function as 

PK{i)^yni{t>0- e.^tK}. 

It is easily seen that 

PK^iS) = sup{(^, p) ; p^K} , 

that is, the gauge function of the polar set coincides with the support function 
of the set K. Let < ci < C2 be such that B^-i(O) C iiT C B^-i(O). Upon 
observing that ^/pk{£.) G K for every ^ 7^ 0, we get 

(2.1) c,\^\<Pk{0<c2\^\: VeeR". 

We say that K S /Cq is of class if dK is of class and all the principal 
curvatures are strictly positive functions on dK. In this case, we define the i-th 
principal radius of curvature at x € dK as the reciprocal of the i-th principal 
curvature of dK at x. We remark that, if K is of class C^, then K^ is also of class 
(see ll^, p. 111]). Moreover, a convex body of class is necessarily a strictly 
convex set. 

Throughout the paper we assume that 

(2.2) K eJC^ is of class C 



2 



Since K will be kept fixed, from now on we shall use the notation p — p^ and 



P° = Pk° ■ 



We collect here some known properties of p and that will be frequently used 
in the sequel. 



Theorem 2.1. Let K satisfy (2.2). Then the following hold: 

(i) The functions p and p'^ are convex, positively 1-homogeneous mR", and of class 
in M" \ {0}. As a consequence, 

pit^^tpiO, Dp{tO = Dp{0, D^pitO^^D^piO, 

p°(U)=tp°(0, Dp"{tO=Dp"{0, D^p\tO^\D'p'{^), 
for every ^ £ R" \ {0} and t > 0. Moreover 

(2.3) {DpiOA) = piO. {Dp'{OA)^P%£,). VeGR"\{0} 

(2.4) i?V(e)e-0, i?V°(O^ = 0, VeeR"\{0}. 

(ii) For every ^,7] £ R", we have 

(2.5) pi^ + v)<PiO+P{v), P'i^ + V)<p"i0+P'iv), 

and equality holds if and only if ^ and r/ belong to the same ray, that is, ^ ^ Xrj or 
rj — for some A > 0. 
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(in) For every ^ ^ 0, Dp{^) belongs to dK°, while Dp°{^) belongs to OK. More 
precisely, Dp[(^) is the unique point of dK'^ such that 

{DpiO, = cind {x, < piO Va: £ /^°, x ^ Dp{£,) . 

Symmetrically, the gradient of Dp'^{^) is the unique point of dK such that 

{Dp°iO, = and {x, < p'iO eK, x^ Dp^'iO ■ 

(iv) The eigenvalues of the second differential D^p at v Q 5""^ are 0, with corre- 
sponding eigenvector v , and the principal radii of curvature of dK'^ at the unique 
point p G dK'^ at which v is attained as an outward normal vector. Symmetrically, 
the eigenvalues of p^ at v ^ S""^^ are 0, with corresponding eigenvector v, and 
the principal radii of curvature of dK at the unique point p S dK at which v is 
attained as an outward normal vector. 

Proof, (i) The convexity and the positive l-homogeneity are a consequence of the 
definition. The regularity in R" \ {0} is proved in H, p. 106]. The identities 
follow upon observing that Dp and Dp'^ are positively 0-homogeneous, whereas D^p 



and D^p^ are positively (— l)-homogeneous. The identity (2.3) is Euler's formula for 



positively 1-homogeneous functions, whereas (2.4) can be obtained differentiating 



(2.3) 



(ii) Follows from convexity of p and p*^ and strict convexity of K and K'^ . 

(iii-iv) See [||, Corollaries 1.7.3 and 2.5.2. □ 

Lemma 2.2. Let K satisfy Then the identities 

(2.6) DpyDp{0)^4^^ Dp{Dp%0) = 



P(0 ' P\i) 
hold for every ^ G M" \ {0}. 



Proof Fixed ^ G M" \ {0}, by Theorem |2]T|(iii) we have that Dp{^) G dK", 
Dp\Dp{i)) G dK, hence p"{Dp{^)) = 1 and 

{Dp"iDp{0),DpiO) = l 

{x, DpiC)) < 1 \fxeK, x^ Dp\Dp{5)) . 



On the other hand by (2^), the point £,/p{£,) G dK satisfies 



which, together with ( p.7\) , implies the first identity in (2.6). The second identity 
can be obtained from the first interchanging the role of K'^ and K. □ 



Remark 2.3. It can be checked that Lemma 2.2 holds under the weaker assumption 
that the sets K, K" G A^g be both strictly convex. This is equivalent to require that 
if is a strictly convex body of class . 

2.4. Nonsmooth analysis. Let us recall the notions of semiconcave function and 
of viscosity solution of Hamilton- Jacobi equations (see for example ^ [l^)- A 
function u: ^4 ^ R is said to be semiconcave if there exists a constant C > such 
that 

\u{x) + (1 - X)u{y) < u{Xx + (1 - X)y) + CA(1 - A)|2; - yf 
for every A G [0, 1] and every pair x,y & A such that the segment [x, y] is contained 
in A. It is easy to check that this amounts to the concavity of the map x i— > 
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u{x) — C|xp on every convex subset of A. As a consequence, if u is semiconcave in 
A then it is locally Lipscliitz continuous in A. 

Let A C K" be an open set, and let w: ^ — > R be a continuous function. The 
superdifFerential D~^u{x) and the subdifFerential D~u{x) oi u dX x ^ A are defined 

by 

„i / s r ™„ , uly) — u(x) — (p, V — x) 1 
D+u{x) = <^ p e M"; limsup — j—^ ^ r ' 

D-uix) = !peR-; hminf "(y) " ^(^) - (P^ ^ - ^) > pi 
[ y^x, j/SA F - 2/1 J 

If u is differentiable at then D'^u{x) = D^u{x) = {Du{x)} (see ^, Lemma 1.8]). 

Let _ff:AxIRxIR"^Rbea continuous function. We say that u £ C {A) is a 
viscosity solution of 

(2.8) H{x,u,Du)^0 in A, 

if, for every x € A, we have that H(x,u{x),p) < for every p 6 D+m(x) and 
H{x,u{x),p) > for every p G D^u{x). 

Example 2.4. Given an open bounded set A C R", the Euclidean distance from 
the boundary of A, defined by 

u(x) = inf \x — y\ x <^ A, 

is the unique viscosity solution of the eikonal equation \Du\ = 1 in yl satisfying the 
boundary condition u = on dA (see p^). 

Let It : A — > R be a locally Lipschitz function. A vector p G R" is a reachable 
gradient of u at x G yl if there exists a sequence {xk)k in yl \ {x} converging to a;, 
such that u is differentiable at x^ for every k and (Du{xk))k converges to p. The 
set of all reachable gradients of u at a; is denoted by D*u{x). 

Since u is locally Lipschitz, it is easily seen that D*u{x) is compact for every 
X & A. Moreover, since u is differentiable almost everywhere in A^ D*u{x) is 
nonempty for every x € A. 

3. Minkowski distance from the boundary 

In the first part of this section we shall define the Minkowski distance from a 
closed set S C R" and we shall prove some basic properties. Then we shall specialize 
to the case S = dil, where C M" is a nonempty bounded open set, and we shall 
study some properties of the Minkowski distance from dfl in Q. 

We recall that if G /Cq is a fixed convex body of class with polar set A'°, 
and that p, p° are respectively the gauge function of K and K'^. 

Throughout this section, S will be a nonempty closed subset of R" . 

Definition 3.1 (Minkowski distance from a set). The Minkowski distance from S 
is the function 5s ■ R" R defined by 

(3.1) Ssix) ^minp^ix -y) , x G R" . 

yes 

Some comments are in order. Let us consider the Euclidean distance from S, 
defined by 

(3.2) d^{x) = mm\x-y\, x G R". 

yes 
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It is clear that dg{x) measures the (Euchdean) length of the shortest segment 
[y,x] joining a; to a point y G S*. In the definition of 6s, the Euclidean distance 
d^{x,y) = \x — y\ is replaced by the Minkowski distance d^^{x,y) = p^{x ~ y). 
Besides the name, the function is not a distance in the usual sense. Namely, it 
satisfies 

d^''(,x,y) > 0, and d^' {x,y) = only li x = y 

d^\x,y)<d'\x,z)+d'^'{z,y) 

for every x,y,z G M", but in general it is not symmetric (unless if is a convex 
body symmetric with respect to the origin). Hence, in the Minkowski distance the 
"length" of a segment [y,x], defined by d^{x,y) — p^{x — y), does not necessarily 
coincide with the length of the segment [x,y]. The space is said to be 

a Minkowski space. It is well known that geodesies in Minkowskian spaces are 
straight lines, so that these spaces are geodesically both forward and backward 
complete (see [|[ Chap. 14]). 

Remark 3.2. As we have underlined in the introduction, we are interested in appli- 
cations to nonlinear PDEs, where the convex body K plays the role of a constraint 
on the gradients of the admissible functions. More precisely, one often deals with 
the set of Lipschitz continuous functions u: A C M" — > IR. such that Du{x) G K 
for almost every x € A, that is p{Du{x)) < 1 for almost every x G A. We shall 
see that the Minkowski distance from S satisfies p{D 6s{x )) = 1 for almost every 
X € A := \ S (see Proposition |3.3| (i), and Theorem 2^(iii)). For this reason we 
prefer to define ^5 in terms of (and hence of K^) instead of using p. 

Since S is nonempty and closed, for every x G K" the set of projections of x in 
S, defined by 

ns{x)^{yeS; 5s{x)^p''{x-y)}, 

is nonempty and compact. 

In the following proposition we gather some basic properties of the Minkowski 
distance 63- 

Proposition 3.3. The Minkowski distance 63 from S is a continuous function in 
R", locally semiconcave in R" \ 5* and, for every x ^ S , the following hold. 

(i) ^5 is differentiable at x if and only iflls{x) = {y} is a singleton; in this 



case DSs{x) — Dp^{x — y) 



anc 



Dp{DSs{x)) = — = — - . 

P{x~y) ds{x) 

(ii) // lis (a;) is not a singleton, then D^5s(x) is the convex hull of the set 
{Dp°{x -y); ye Usix)}, while D'Ssix) = 0. 

(iii) For any y G 115(2;), and for any point z in the segment {x,y) we have 
115(0) = {y}- Hence ds is differentiable at every point z G {x,y) and 
D5s{z) = Dp\x-y). 

(iv) D*5s{x) = {Dp"ix ^y); ye Us{x)}. 

Proof. The proof below is a straightforward adaptation of Proposition 2.14 in 
and Corollary 3.4.5 in where the case K = Bi{0) is considered. 
The continuity of 6s follows from the inequalities 

-P°{z ~x)< 6s{x) - 6s{z) < p°{x - z), Vx, z, G R" . 
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The other properties follow from the general theory concerning the marginal func- 
tions (see H Proposition 2.13], ||, Proposition 3.4.4]). Let A = M" \ 5 and 
F{y,x) — fr{x — y), so that 6s{x) = min^gg F(?/, x) for every x £ A. Since p° 
is of class in M" \ {0}, it is clear that F{y,x), DxF{y,x) = Dp°{x — y) and 
DxxF(y, x) = D^p'^{x — y) are continuous functions in S" x A. The local semicon- 
cavity of Ss in M" \ S" is a consequence of the same property for marginal functions 
proved in [|| Proposition 3.4.1]. 
If we define the sets 

M{x) := {y e S; Ss{x) = F{y, x)} = \ls{x) , 

Yix) := {DxF{y, x); y G M{x)} = {Dp\x ~ y); y e Ugix)} , 

propert ies ( i) and (ii) are an immediate consequence of Proposition 3.4.4] and 
Lemma 2.2 . 

In order to prove (iii), let us start by proving that y G 115(2) for every z G {x, y). 
Namely, if not, there exists y' & S such that p^[z — y) > p^{z — y'). Since z G {x, y), 
by Theorem |]^(ii) we obtain 

p'ix - y) = pO(x -z) + p°{z -y)> p\x - z) + p°{z - y') > p\x ~ y') , 

which contradicts the fact that y G fig (a;). The same argument shows that ^s{z) ~ 
{y}. Namely, if we suppose that there exists y' G \ls{z), y' ^ y, we have that x — z 
and z ~ y' are not proportional, hence, again from Theorem |2.lK ii), we get 

pO(x - y') < p°(x -z) + p\z ~ y') = p\x - z) + p\z - y) = p°(x - y), 

which contradicts y G Iis{x)- Now the differentiability of 83 at z and the fact that 
D5s{z) = Dp^{x — y) follow from (i) and from the 0-homogeneity of Dp^ . 

Finally, in order to prove (iv), we have to show that D*Ss{x) — Y{x). The 
inclusion D*Ss{x) C Y{x) follows from Proposition 3.4.4]. Now let p G Y{x), 
and let us prove that p G D*dnix). By definition of Y{x), there exists y G Ils{x) 
such that p = Dp'^{x — y). From (iii) it we have that 5s is differentiable at any 
point z G (x, y), and 5s{z) = p, hence p is a reachable gradient at x. □ 

The following lemma, which states that the multifunction x 1— > T\s{x) is sequen- 
tially upper semicontinuous, is a straightforward consequence of the continuity of 
p° and 5s- 

Lemma 3.4. // {xk)k C M" is a sequence converging to a point x G M", if yk G 
T\s{xk) for every k, and if the sequence {yk)k converges to y, then y G 1X5(2;). 

Proposition 3.5. Let x ^ S . For y G Iis{x) let us define 

(3.3) ly{x) = Ss{x) ■ sup{t >0; ye Hsiy + t{x ~ y))} G (0, +00] , 
and, if ly{x) < +00, the point 

(3.4) m,y{x) = y + ly{x) ^^—^ . 

5s{x) 

Then ly{x) and my{x) (when ly{x) is finite) do not depend on the choice of y e 

Proof. Observe that, if z = y + t{x ~ y), t > 0, and y G Hs{z), then 
5s{z) = p\z - y) ^tp°{x - y) ^t5s{x) . 
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From the very definition of ly{x) and Proposition 3.3(iii), we have that ly{x) > 
Ss{x), for every y £ ns'(a;). We have two possibihties: 

(i) ly{x) = Ssix) for every y G Ils{x); 

(ii) ly{x) > Ss{x) for some y G Ils{x). 

In case (i) it is plain that ly{x) does not depend on the choice of y G 115(0;) 



Furthermore, the right hand side of (3.4) coincides with x for every y G n5'(a;), 
hence my{x) = x is also independent of the choice of y G Ils{x). Let us consider 
case (ii). Let y G Ils{x) such that ly{x) > Ss{x). By definition of ly{x), there 
exists t > 1 such that, setting z — y + t{x — y), y & Ils{z). Since x G (y, z), from 



Proposition 3.3(iii) we have that Hsix) ~ {y} is a singleton. This concludes the 



proof of the proposition. □ 
Thanks to Proposition |3.5| , the function / : M" \ 5 ^ (0, +00], 

(3.5) lix) 6six) ■ snp{t > 0; y G Hsiy + t{x - y))] (y G IVs{x)) 

is well-defined. The same holds true for the function m: {x G R" \ 5; l{x) < 
+00} ^ defined by 

(3.6) mix):^y + l{x)^^ {yeUsix)). 

Remark 3.6. It is straightforward to see that, given x ^ S and y G 115(0;), 

(3.7) r(a;)=sup{A>0; yGn5(a;A)} , xx.= y + \ 



pO{x - y) ■ 

It is easily checked that x\ = x for A = Ss{x), hence, from Proposition |3.3| (iii), 
ns(a;A) = {y} for every A G [0, 6s (x)). As a consequence, we recover the inequality 
i{x) > 6s{x). 

From now on we shall restrict our attention to the case S — R" \ 51, where 
51 C R" will be a fixed nonempty connected bounded open set. 

Definition 3.7 (Minkowski distance from the boundary). The Minkowski distance 
from the boundary of Q is defined by 

(3.8) do(a;) = ^R,>\o(a;) = inf p°{x-y), a; G . 

yedU 

For every a; G fl we shall denote by Tl{x) ~ 119^(2;) the set of projections of x in 
951. For every xq G 951 we define the proximal normal cone to 951 at a;o as the set 

(3.9) N^{xo) := {Xv; xq + v e 51, xq G n(a;o + v), A > 0} 
(see ll^). In general, N^{xq) is convex but it may not be closed. 

Definition 3.8 (Singular set). We say that a; G 51 is a regular point of 51 if Il{x) 
is a singleton. We say that a; G 51 is a singular point of if a: G 51 is not a regular 
point. We denote by S C 51 the set of all singular points of fi. 



From Proposition 3.3(i), S coincides with the set of points in 51 at which dn is 
not differentiable. 

Since 51 does not contain half spaces, it is clear that, for every x E CI, the quantity 



l{x) introduced in ()^^ is finite, hence the point m{x) (see (|3.6|)) is well defined. 

Given x € fl, we call a ray through x any closed segment [y, m{x)] where y G 
n(a;). Thus there is a unique ray through a; G 51 if and only if n(a;) is a singleton. 
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that is, if and only if a; S. From Proposition 3.5 it is clear that, in any case, all 
rays through x £ H. have the same length l{x) in the Minkowski distance. Observe 
that the length of a ray [y, rn(a:)] is measured by p^{rn{x) — y) and, in general, it is 
different from the length p'^(y — m{x)) of [m{x), y\. 



In Lemma 3.9 below we prove that / is an upper semicontinuous function in VL. 



Furthermore, in Proposition 3.1C we shall define an upper semicontinuous extension 
Z of Z to Q,. 



Lemma 3.9. The Junction ^ M defined in (3.5) is upper semicontinuous in 

n. 

Proof. We follow the lines of the proof of Proposition 3.7 in [|l3j, where the same 
property corresponding to the Euclidean distance is proved. By contradiction, 
assume that there exists a sequence {xk)k C il converging to a point x G such 
that 

(3.10) limf(xfe) > [(a;). 

k 

For every A; G N, choose yk G Il{xk) and define 

, . [ iM „ lix) ] 

Afc = mm < -— — - , 2- 



dn{xk) ' dn{x) j 

The sequence {yk)k C dil is clearly bounded. Moreover, since l{xk) > dn{xk) 



for every fc G N (see Remark 3.6), we have 1 < Afe < 2 l{x) / dn{x) , hence the 
sequence {Xk)k is also bounded. We can then extract a subsequence, which we do 
not relabel, such that limfe Afc = A > 1, linife yfe = y E dfl. Since, by definition of 
l{xk), Afc < l{xk)/dn{xk), we have that yk € Tl{yk + ^ki^k - yk)) for every A: e N. 



Hence, by Lemma 3.4, we infer that y G n(y + A(a; — y)). By the definition of l{x) 
we conclude that A < l{x)/da{x), so that A^ = l{xk) / du{xk) for k large enough. 
Then 

, l{xk) 

k dfi(Xk) 

that is 

liml{xk) = A \imdn{xk) < lix), 

k k 



in contradiction with ( 3.1C| ). □ 



Proposition 3.10. The function I: ft ^M. defined by 

J l{x) if X G il, 

]s\rp{dfi{z); z Eft and x G n(^)} if x G 9f2, 



(3.11) l{x) 



is upper semicontinuous in Q. Here I: — > (0, +oo) is the function defined in (S.L) 



Proof. Recalling Lemma 3.9 we have only to prove that I is upper semicontinuous 
at every point of dQ. Let xq G dfl. Assume by contradiction that there exists a 
sequence {xk)k C fl, converging to xq, such that 

(3.12) \iml{xk) > Kxo) ■ 

k 

Upon passing to a subsequence, it is enough to consider the following two cases: 

(a) Xk G dfl for every fc G N; 

(b) Xfe G for every A; G N. 
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In case (a), by the very definition of I, for every A; e N there exists a point Zk such 
that 

(3.13) ZfcGn, Xfcen(zfc), dn{zk)>l{xk)-l/k (fc > 1) . 

Since {zk)k is a bounded sequence, there exists a subsequence, which we do not 
relabel, and a point Zq G ^ such that linife Zk — zq. From the upper semicontinuity 



of n we have that xq £ n(zo), hence by ( |3.13| ) 



l{xo) > dn{zo) > limsupZ(a;fe) , 

k 

in contradiction with ( |3.12| ). 

Consider now case (b). For every fc G N, let yk G Il{xk)- Up to a subsequence 
we can assume that {yk)k is convergent. By the upper semicontinuity of 11, {yk)k 
converges to xq, since n(a;o) — {xq}. Let us define 

a= or''""' V 

p"(xk - Uk) 



By the characterization (3.7) of I, for every fc we can choose Afc > l{xk) — 1/fc, 
\k > 0, such that the point Zk — yk + ^k£.k satisfies yk G Il{zk)- Since {Xk)k is 
bounded and S^k G for every fc, we can extract another subsequence (which 

we do not relabel) such that lim^ ^fc = ^ G dK^ and limfe = A. Then we have 
limfe Zk = Xq + and, by the upper semicontinuity of 11, xq G TI{xq + A^), hence 
K^o) > dn{xo + A^). Collecting all the information we obtain 

li^o) > dn{xo + XS,) = X = limAfe > limsup;(xfe). 



which contradicts ( 3.12| ). □ 



Definition 3.11. The point m(x), defined in (3^), is called the ridge point of 
X G il. The set 

TZ = {x eTl; dn{x) = l{x)} 

is called the ridge of Q,. 

Observe that, if x G 7?., then either x — m{x), or a; G dfl and x ^ n(?/) for 
every y G 51. If a; G il is not a ridge point, then l{x) > dn{x), and, from Proposi- 



tion 3.2 (iii), the set n(a;) is a singleton, hence dji is differentiable at x. Thus the 
singular set E is contained in the ridge TZ of fl. We remark that, in general, the 
two sets do not coincide, as it is shown in the following example. 

Example 3.12. Let fl = {{x,y G x^/a'^ + y^/P < 1}, where < 5 < a and 
K = Bi{Q). The points P = ((a^ - b'^)/a,0) and Q = (-(a^ - b'^)/a,0) are the 
centers of curvature of dfl at (a, 0) and (—a, 0) respectively. It can be checked that 
n = [Q, P] whereas S = (Q, P). 

4. Distance from the boundary of a smooth set 
Throughout the rest of the paper, we assume that 
(4.1) n C K" is a nonempty, bounded, open connected set of class C^. 

For every xq G dfl we denote by ki{xo), . . . , k„_i(xo) the principal curvatures 
of dQ at Xq, and by i^^xq) the inward normal unit vector to dQ at xq- 

Since the boundary of fl is regular, we will be able to extend dn outside fl in 
such a way that this extension turns out to be of class in a tubular neighborhood 
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of dft (see Theorem i.ld below). This fact will allow us to define Ddn and D^dn 
on points of dH,. The extension of dn to R" can be constructed as follows. Let us 
define 

(4.2) d^{x)= inf pO(y-x), xeW^\n, 
and the signed distance function from dfl 

(4.3) dUx)^h^^}\ ifxen 

For every x £ M" \ define the set of projections from outside fl 

n-{x)^{yedn; d^{x)^p"{y-x)}, 
and we extend the projection operator 11 to be 

^Tl{x), ifxGn, 



U^ix) = 



n"(.T), ifxeM"\f7. 



We define also the set of singular points of d^ in M" \ fl and the set of 
singular points of d^^ ™ I^"- We clearly have the inclusion C E U E^ U dfl. 

It is worth to observe that, in general, 11" (x) 11^^(2;) for x G M" \ fi, since /9° 
need not be symmetric. 

It is clear from the definition that d^ is the Minkowski distance from dfl induced 
by the gauge function of —K^. We can interpret dQ{x) as the Minkowski distance 
of a point x € R" \ to the boundary of fl. 

As in ( 3.11 ), we define a function on dfl by 

(4.4) r{xQ)^sup{d^{z); zeMJ'X n and xoeIl~{z)} xq e dn . 

In the following lemma we show that, under our regularity assumption on 17, the 
functions I and are bounded from below by a positive constant on dQ. 



Lemma 4.1. Let Q satisfy (4.1). Then there exists a positive constant ji such that 
1{xq), {xq) > fi for every xq G dfl. As a consequence, TZ Cfl. 

Proof. We shall only prove the assertion concerning I, the other being similar. 
Since is of class , then it satisfies a uniform interior sphere condition of radius 
r > 0, that is, for every boundary point xq G dVl there exists zq G 17 such that 
Xq G Br{zo) C il. Let i? > denote the maximum of the principal radii of 
curvature of so that if" (and —K'^) slides freely inside Br, that is, for every 
boundary point y G dBn there exists z G B^ such that y E z — K'^ C Bn (see 



23, Corollary 3.2.10]). We are going to prove that, for /i = r/R, —pK^ slides 



freely inside VL. Let x^ G dO,. We have to show that there exists z G 57 such 
that Xq a z — ji C f7. From the uniform interior sphere condition, there exists 
zq E ^ such that G i3r(^o) C 17. On the other hand, — /i A'" slides freely inside 
hBr = Br- Thus, there exists z G Br(zo) such that xo G z — /i AT" C Br{zo). We 
have proved that for every xq G 917 there exists z G 17 such that G z — /iAT^ C 17. 
Then p = p°(z — xo) < p^{z~y) for every y G 917, that is xq G n(z) and (ia(z) = /i. 
Hence 1{xq) > p.. □ 



Remark 4.2. As a consequence of Lemma 4.1, we obtain that, for xq G 917, the 
supremum in the definition ( |3.11 ) of l{xa) is achieved by a point z G 17. Namely, 
from the continuity of dn and the upper semicontinuity of 11, the supremum is 
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achieved at a point 2 G il. Should z £ dil, then xq G n(z) would imply z — xq and 



l{xo) — 0, in contradiction with Lemma 4.1 



In the definition (4.4) of I (xq) the supremum need not be achieved, but the 
argument above shows that we can replace M" \ by M" \ 51 in the definition of 

The following result relates the gradient of to the inward normal of dft. 

Lemma 4.3. Let a; G M" \ E^ a; ^ dil, and let n"(a;) = {xq}. Then 

Proof. Consider first the case x G f2 \ S. From Lemma 3.5 in we have that 
Ddit{x) = Xi'{xq) for some A > 0. On the other hand, fro m Proposition 3^(i) we 
have that Dd^{x) = Dp^{x — xq), and, from Theorem |2.l| (iii), this vector belongs 
to dK. Then 1 = p{Ddi-i{x)) = \p{i^{xo)), that is A = l/p{iy{xo)) and follows. 
Consider now the case x G K" \ (SI U From Proposition ^^{i) and Lemma 3.5 
in we have that _Dp(_^-)o(x — xq) = Dd^{x) — —Xv^xq), for some A > 0. Upon 
observing that 

we deduce that Ddf^{x) — Dp^{xQ — G dK , hence 

1 = p{Ddt,{x))^\p{iy{x^)), 
and (|^ follows. □ 

We are now in a position to characterize the normal directions to the boundary 
oin. 

Proposition 4.4. For every xq G dfl, the proximal normal cone ofVl at xq, defined 
in ^3.d( ), is given by 

(4.6) 7V^(xo) = {A Dpiiyixo)); A > 0} , 

whereas the proximal normal cone o/K" \ il at xq with respect to is —N^{xo). 



Proof. By Proposition 3.3(i) and Lemma 4.3 the vector Dp{i'{xo)) belongs to 
N^{xo). On the other hand, by the definition of N-^{xo), we have that w G 
Af^(a;o)\{0} if and only if there exists ^ > such that xq G II{xo + pw). Then, for 
e G (0,^), from Proposition |3.3| (iii) we have n(a;o + ew) = {xq}. From Lemma 4.3 
and Proposition |]3|(i) we have that 

w 

Dp{v{xq)) = Dp{Ddn{xQ + „, , , 

p"('W) 

hence w = \Dp{v{xq)) with A = p^(w) > 0. The computation of the proximal 
normal cone of M" \ at is similar. □ 



Remark 4.5. As a consequence of Proposition 4.4 and Proposition |3.3| (i) it is clear 
that, from any point Xq G Dp(y{xQ)) is the unique inward "normal" direction 
with the properties 

li^lxQ + tDp(z/(xo))] = {xq}, df^ixo +tDp{iy{xo))) = t 

for t small enough (see also Lemma 2.2]). In the following proposition we shall 
show that, under our assumptions, the relations above hold for t G (— ^^(xo), 1{xq)). 



This is a well known fact in Riemannian geometry (see |22, §111.4 
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Proposition 4.6. For every xq G dfl, we have that 

(4.7) n'{xa+tDp{iyixo)))^{xo} e [-1' ixo),l{xo)), 

whereas xq ^ n^(a;o + t Dp[v{x{)))) for every t < — /^(xq) and for every t > 1{xq). 

Proof. We prove the assertion for t > 0, the case t <0 being similar. Fixed xq G dil, 



Proposition 



by Remark 4.2 there exists z G such that xq € 11(2:) and d^^ixo) = 1{xq). By 
Proposition 



3.3(iii), {xq} — H{x) for every x in the segment {xq,z). Moreover, by 



3.3 (i) and ([4.5|), the segment [xq, z) can be parameterized by 



Xt = Xq +t-^ = +tDp{v{xQ)), t G {O,l{xo)) , 

p"{z - Xo) 



SO that (4.7) holds true. On the other hand, by the very definition of 1{xq), it 
cannot happen that xq G 11^(2:4) and t — dn{xt) > l{xo)- O 

The following result, which will be used in the sequel, states that every point 
a; G fi \ S belongs to the interior of the ray {p{x), m(x)), where p{x) is the unique 
projection of x on dH,, and m{x) is the ridge point of x defined in (|3.6|). 



Lemma 4.7. Let a; G 51 \ S and let p{x) denote its unique projection on dfl. Then 
there exists a point xi G 51 \ S such that n(a;i) = {p(x)} and x G {p{x),xi). 

Proof. Since SI \ E is an open set, there exists r > such that Br{x) C SI \ S. For 



every z G Br{x) let p{z) denote the unique projection of z on dfl. From Lemma 3.4 



the map p is continuous in Br{x). Let V C dVL be a local chart on 9S7 containing 
p{x), and let 1^: y — > R"^^ be a local coordinate system, that is, a bijection 
from V to the open set U — ip{V) C M"^^. Since p is continuous in Br{x), there 
exists 5 G (0,r) such that p{Bs{x)) C V. Consider now the map V": 5*""^ R"^^ 
defined by 

iiiu) = •f{p{x + 5 u)) , ugS"""^. 

It is clear that this map is continuous, being the composition of continuous maps. 
From the Borsuk-Ulam theorem (see [|ll|. Corollary 4.2]) there exists uq G 5*"^^ 
such that ip{uo) = ip{—uo). Since ip is one-to-one, we deduce that p{x + Suq) — 
p{x — Suq), that is, the points xi ~ x + Suq and X2 = x — Suq have the same 
projection xq G From Proposition |3.3| (i) and Lemma L3 we have that 

Xi = Xq + dnixi)Dp{i/{xo)) , i = l,2. 

It is clear that the points xq, x\ and xi are collinear. Assume, just to fix the 
ideas, that d^(x\) > dn{x2), that is, p°{xi — Xq) > p°{x2 — Xq), in such a way 
that X2 G {xo,xi). Since x G (0:1,0:2) by definition of xi and X2, we also have 
X G {xq,xi). Finally, recalling that n(a;i) = {xq}, from Proposition p.3K iii) we 
conclude that p{x) — xq. □ 



Remark 4.8. As consequence of Lemma 4.7, we have that l{x) > dn{x) for every 
a; G fl \ S. Hence 7^ C S. 

For many calculations wc shall use a preferred system of coordinates in order to 
parameterize dVt in a neighborhood of a point xq. 

Definition 4.9 (Principal coordinate system). Let xo G dVl. We call principal 
coordinate system at xq the coordinate system such that xq — 0, e„ = iy{xQ) and 
Ci coincides with the ?-th principal direction of Oil at xq, i = 1, . . . ,n — 1. 
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Using the principal coordinate system at xq, dfl can be parameterized in a 
neighborhood of xq by a map 

where U C W^~^ is a neighborhood of the origin, and (f): U ^ M is a map of class 
satisfying 

(4.8) 0(0) = 0, 1^(0) -0, i=l,...,n-l. 

We shall refer to X as a standard parametrization of dfl in a neighborhood of xq. 
If we denote by N{y) — v{X{y)), y then 

Vi + \Dy(f>{y)\- 

where Dy<j) = (d(f)/dyi, . . . , dcjj/dyn-i). Since, in the principal coordinate system, 
the vectors ei, . . . , e„_i are the principal directions of dfl at xq, we have 

dN 

(4.9) iV(0) = e„, — (0) = -K,e„ z = 1, . . . , n - 1, 

oy, 

where ki,...,k„-i are the principal curvatures of dfl at xq. Starting from the 
identities 

^^=-N,{y)^l + \Dy(t>{y)\^, j = 

and differentiating, we easily obtain 

dN- 

(4.10) — 1-(0) = -— ^(0) = «,,5,„ = 

oyioyj oyi 

where 5ij denotes the Kroneker symbol. Moreover, since 

for every j, j = 1, . . . , n — 1, by (4.8) and ( 4.1Cl| ), we get 

dX d^X 
(4-11) -^(0) = ei, — — — (0) = Ki(5y e„, i, j = 1, . . . ,n - 1 . 

oyt oyioyj 

In the following lemma we introduce the main tool of our theory, that is a 
parametrization of that will allow us to prove the regularity of the signed distance 
near dfl, as well as the regularity of the set S. Moreover, in Section]^ we shall use 
this parametrization in order to prove a change of variables formula in multiple 
integrals on fi. 

In the following computations, we extend the notion of inward normal for x in 
a tubular neighborhood of dfl, by setting iy(x) to be the gradient of the Euclidean 
signed distance of x from dfl. In this way Dv{xo) is well defined for every xq £ d^, 
and Dv{xf))v(xo) = 0. 

Lemma 4.10. Let Y : U ^ E", U C K"~^ open, be a local parametrization of dfl 
of class . Let U xR be the map defined by 

(4.12) ^{y,t)=Yiy)+tDpiiy{Yiy))), iy,t)eUxR. 

Then ^ E C^{U xM), and 



(4.13) det D^{y, t) = pHY{y)) yW) det[/ + t D^p{v{Y{y))) Dv{Y{y))] 
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for every {y, t) ^ U x where 

I dY dY \ 

g{y) ^ det{gij{y)), gij{y) = {^ — {y), —{y) j , i, j = 1, . . . , n - 1, 
is the determinant of the matrix of the metric coefficients. 

Proof. Since Dp e (M" \ {0}) and o F e C^W) , 4' is of class in U xR. Let 
us fix yQ e U, let = Y{yQ), and define 

dY 

Q = I + tD'^p{v{xo))Dv{xo), = — (2/0), i = 
We have that 

detD*(?/o,0 = det [Q wi, . . . , Q w„_i , L»p(i^(a;o))] ■ 

Notice that, for every i = 1, . . . , n — 1, {wi, v{xq)) = and, thanks to ( ^.4| ), 

(L'^p(i/(a;o))£'j^(a;o)t«i, t^(a;o)) = {Diy{xo)wi, D"^ p{iy{xQ)) v{xo)) =0. 

Hence {Qwi, iy{xQ)) = for every i — 1, ... ,n — 1, and Q iy{xo) — iy{xo). Upon 
observing that Dp{v{xq)) = w + p{v(xq)) v{xq), with {w, v{xq)) — 0, we get 

<ieiD'^{yQ,t) ^ p{v{xo)) det wi, . . . , Q w„_i , i^(xo)] 

= Pii'ixo)) det Q det [wi, . . . , Wn-i,v{xo)\ . 

Finally, the matrix C — [wi, . . . , ^{xq)] satisfies 

(C^ C)y = (wi, Wj) , i,j = l,...,n- I, 
(C^ C)fe„ = (C^ C)„fc = 4«, fc = 1, . . . , n . 

Hence (det C)^ = det(C'^ C) = ^(yo), and is proved. □ 

The crucial point in what follows is to have exact information about the degen- 
eracy of the map 5". The first step in this direction is to simplify the computation 
of det[I + tD^p{iy)Du]. 

Lemma 4.11. For every xq £ dfl, we have that 

(4.14) det[I + t piiyixo)) Di^ixo)] = det(/„_i -tRD), 
where D, R E A4n-i are the symmetric square matrices with entries 

(4.15) = -Ki(a;o) % , Rij ^ {d'^ p{v{xo)) ej, ei) , i, j = 1, . . . , n - 1, 
in the principal coordinate system at xo . 

Proof. It is enough to observe that, in the principal coordinate system at xq, 
{Du{xq) Ci, Cj) — —Ki 6ij for every i, j — I, . . . ,n ^ I, whereas, from 

(D"^ p{iy{xQ)) Diy{xo) e^, j^(xo)) = {Diy{xo) e^, D^p{v{xo)) ^ixo)) = 

for every i = l,...,n — 1, hence ( [1.14[ ) follows. □ 



The equality (4.14) is fruitful, since it allows us to deal with a positive definite 
matrix R with known inverse matrix, as stated in the following result. 



18 



G. CRASTA AND A. MALUSA 



Lemma 4.12. Let xq G dfl, and let H G Mn-i be the matrix with entries 



(4.16) H,, - (D^p'^iDpii^ixo))) ej, e,) , i,j 
in the principal coordinate system at xq . Then 

(4.17) p{v{xq))RH ^In^^, 



1, 



1, 



where R G A^n-i is the matrix defined in (^.15). As a consequence, R and H are 
positive definite matrices. 



Proof. In order to prove (4.17), let us differentiate the first identity in (2.6) with 
respect to ^, obtaining 

1 ^ 1 



I - 



■ c^Dpio, ve^o. 



From the symmetry of D^p'^ and D^p, the adjoint matrix identity can be written 

1 1 



as 



D'p{OD'p"{Dp{0) = 



I - 



■DpiO^^, Ve^O. 



P(0 ^ piO' 

Notice that, if v is orthogonal to = i^{xo), then {Dp{v) O z^) u = 0, hence, in the 
principal coordinate system, for every i, j = 1, . . . , rt — 1, 

^ ■5,, = {D^p{v)D^p\Dp{v))e,,e,) 



P{^) 



{D'pH eu e,) {D^p'^iDpiv)) e„ ei) , 



1=1 



where the last equality is due to the fact that, from (2^), p{v) e„ — D^p{v) v — Q. 
The relation above is exactly ( 4.17 ). 

Recalling that D^p(y) and D^p^{Dp{v)) are positive semidefinite (due to the 
convexity of p and p ^), it follows that R and H are both positive semidefinite. 
As a consequence of ( 4.17 ), R and H are also invertible, hence they are positive 
definite. □ 

The fundamental step in the proof of the regularity of in a tubular neigh- 
borhood of dVl and of do in \ S is the following. 

Theorem 4.13. Let Y : U ^ M", U C M"~^ open, be a local parametrization of dfl 
ofclassC^. Let^: UxM. R" be the map defined in Then det D"^ {y , t) > 

for every y & U and every t G {—l~{Y{y)),l{Y{y))). 



Proof. Let us fix xq G 30.. There is no loss of generality in assuming xq = Y{0). 
In the principal coordinate system at xq, we have also = 0. We shall denote 
V = z^(xo) = e„. From Lemmas 4.10 and 4.11, in this coordinate system we have 
to prove that 

dct{In-i-tRD) > 0, Vi G {-1-{xq),1{xo)), 



(4.18) 

where D, R E AAn-i are the square matrices defined in (4.15). 
Observe that, from (4.17) 



Ln-i-tRD = p{v)RH -tRD ^ p{v)R [H- 



D 



where H is the matrix defined in ( 4.16| ). Since det(i?) > by Lemma 4.12, in order 
to prove (4.15), it is enough to prove the following claim. 
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Claim. The matrix A{t) G Mn-i defined by 

A{t)=H ^D, <gM, 

is positive definite for every t G {—1^{xo),1{xq)). 

Proof. For t = we liave ^(0) = H, which is positive definite by Lemma 4.12, We 
shall now prove the claim for t G (0, ^(a;o)), the case t G {—l{xo), 0) being similar. 
We start proving that, for every t G {O,l{xo)), A{t) is positive semidefinite. Let 
zt = xq + tDp{v). From Proposition 4.6 we have that n(zt) — {xq}, and dn(-Zi) = 
t = p^{zt — Xq). Moreover, the following ball in the Minkowski norm 

B = {xe M"; p°{zt~x) < t} 

is contained in il, and xq G dB n dfl. Thus 

h{y):=p'>izt-X{y))>t Vj/ G h{0) = t . 

Then, the function h G C^{U) has a local minimum point at y = 0, hence its Hessian 
matrix at y = must be positive semidefinite. It is straightforward to check that 

— — (y) = ( D'p'^izt - — (y), — (y) 

oyioyj \ oyj oyi 

(0) = {D'^p°{zt - Xq) e^, Cj) - Ki Sij {Dp°{zt - xo), v) • 



Using (4.11) we obtain 
dyidyj 

Since zt — xq — tDp{h'), from the positive 0- homogeneity of Dp'^, the positive 

2„0 



(— l)-homogeneity of Dp, and ( |2.6| ) we get 

Dp"{zt - xo) = Dp^{Dp{uj) = - 2:0) = y D''p^{Dp{u)) , 

Piny)) t 



hence 



-iO) {D'pOiDpi.))e., e,) - ^S, 



that is, D'^h{0) = {l/t)H - {l/p{v)) D = A(t)/t. Since the Hessian of /i at y = is 
a positive semidefinite matrix, we conclude that A(t) is also a positive semidefinite 
matrix for t G (0, ^(xo)). 

Let us prove that, in fact, A{t) is positive definite for every t G (0, ^(a;o)). It 
is plain that it is equivalent to prove that A piv) H — D is positive definite for 
every A > l/l{xo). Let us fix A > l/l{xo) and choose A' G (l/^(a;o). A), so that 
A' p(z^) H ~ D is positive semidefinite. Since, by Lemma ^.12| , H is positive definite, 

((A H - D)w, w) = ((A' p(i/) H - D)w, w) + (A - A') H w, w) 

> {\-\'){p{v)Hw, w) > 

for every w 7^ 0, that is, A p{v) H — D is positive definite. This concludes the proof 
of the claim and of the proposition. □ 
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Remark 4.14. It is well known that is the unique viscosity solution of the 
Hamilton- Jacobi e quat ion p{Ddn) = 1 in vanishing on dfl (see Q ||, |l^). The 
map ^' defined in (4.12) gives the characteristic curves associated to this PDE. The 
regularity of the map proved in Theorem 4.13 will be used in Theorem 4.16 in 



order to prove the regularity of df^ in a tubular neighborhood of dil. It should 
be noted that this kind of result in fact follows from the local existence theory for 
first order PDEs based on the method of characteristics. Nevertheless the regularity 
of the map 5" is an essential tool in order to prove the regularity of dji in the 
whole set \ E (see Theorem 3.10 below). 

Corollary 4.15. The ridge set TZ has zero Lebesgue measure. 

Proof. Let Yk'.Uk ^ K", Uk C R"^^ open, k ~ 1, . . . , iV, be local parameterizations 
of dVl of class C^, such that [_}k=iYk{Uk) = dVL. For every k = l,...,iV, let 

*fc : Wfe X K ^ M" be the map 

^k{y,t) Yu{y) + tDp{u{Yk{y))). {y,t) eUk x M. 

For every fc = 1, . . . , n let [7^ C Z//fc be a compact set such that IJ^, Yk{Uk) covers 
dft, and let 

Ak^{{y,t); yeUk, t e [OJiYkiy))]} . 
From Proposition 3.1C, / is an upper semicontinuous function, hence for every k — 
1, . . . , n, Ah is a compact set and the Lebesgue measure of the graph 

^^\n)nAk - {iy,t) eUkXR; t = liYkiy))} 

vanishes. From Theorem [4 . 1 3| we know that, for every k — I, . . . , N, ^ C^{Uk x 
K), hence it is Lipschitz continuous on the compact set Ak. Let L be the maximum 
of the Lipschitz constants of the functions ^'i, . . . , ^at. Since ljfc=i ^fc(^fc) = 
and hence 7^ C Ijf^^ {^"^{TZ) n ^fe), we finally get 

N N 

< J2 ^" ["^k '(7^) n a^)] < ^ [*fc '(7^) n a^] = o 
fe=i fc=i 

and the proof is complete. □ 

Theorem 4.16 (Regularity of df^). The function df^ is of class in a tubular 
neighborhood of dil of the form 

Af,^ {x e M"; -fj. < d^{x) < p,} 

for some fi > 0. Furthermore, for every xq G dQ 

(4.19) Ddf,{xo + tDp{:y{xo)))^Ddf,ixo), te{-p,p), 

and then the identity (^.5) holds also for x £ // O is of class C'^'" and p 6 

C'='"(]R" \ {0}), for some k > 2 and a € [0, 1], then df^ is of class C"='" in A^. 



Proof. From Lemma 4.1 there exists p, > such that l{xo),l (xq) > p for every 
Xq € dfl, hence, thanks to ( |4.7| ), for every point z G A^ the projection Il^{z) is a 



singleton. Let zq G A^j^ and prove that df^ is of class in a neighborhood of zq. 
Let n(zo) = {^o}, and let Y : U R", U C K."^^ open, a local parametrization 
of dft in a neighborhood of xq with Y{0) = xq. Let E C^{U x R) be the map 
defined in ( [4.12[ ), and let to = ^^(20) G (~M)M)- From Theorem 4.13 we have that 
det £'\E'(0, to) > 0, hence from the inverse mapping theorem it follows that there 
exists a neighborhood V C of zq such that the maps y — y{x) and t — t[x) are 
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of class C^{V). Since dfJx) — t{x) for every x (1 V, this proves that df^ is of class 
C^{V). Moreover, by Ddf^{^{0,t)) = v{xo)/ p{iy{xo)) for every i £ (0,to), 

hence by the continuity of Ddf^ in and the fact that (4.1£) holds in {0,to), we 
have 

(4.20) Dd'^{xo)^Hxo)/piHxo)) (xoedQ). 

Finally, from ( |4.5| ) we have that Ddn{x) ~ N{y{x))/ p{N{y{x))), x G V, where 



iV = jy o r e C\U), hence Dda G Ci(T/), that is dn £ C2(F). 
The last part of the proposition follows from the fact that, if fl is of class C*^'" and 
p e C'''"{R" \ {0}), then * e C^-^^°'{U x M). □ 



Remark 4.17. As a consequence of Theorem 4.16 , we have that C R" \ dVt. In 
particular E C and C M" \ H. 

We conclude this section giving an explicit representation of the Hessian matrix 
of dfj in xq G do. with respect to the principal coordinate system. From now on, 
we shall denote D'^dQ{xo) ~ D^df^{xo) for every xq S dfl. 

Lemma 4.18. Let xq G dCl. Then, in the principal coordinate system at xq we 
have that 

1 A ■ ■ 



(4.21) 



(4.22) {D^daixo)e„ iy{xo)) 

(4.23) (D^dnixoHxo), i^{xa)) 



1, 



, n — 1. 



{Dp{iy{xo)), e^) , i = 1, 

n-1 



, n — 1. 



1 



^ K,, {Dp{iy{xQ)), CiY 



(Recall that i'{xo) — e„ m the principal coordinate system.) 

Proof. If we differentiate the identity dQ{X{y)) — for every y & U with respect 
to t/i, i = 1, . . . , n — 1, we obtain 



r)X 

Ddn{X{y)\ —iy) 



0, yeU. 



A further differentiation with respect to yj, j ~ 1, ... ,n — 1, gives 

D^d.iXiy))-iy), -iy)) + (l...(X(,)), ^(,) 

Ev alua ting the last expression at y = 0, and taking into account ( 4.11 ) and Lem- 
ma O, we obtain ( 4.21 ) 



0, yeU. 



In order to prove ( [4. 221 ) and ( ^4.23D , let us differentiate the identity ( ^4.191 ) with 
respect to t, and evaluate the result at t — Q. We obtain 

(4.24) D^dnixo)Dpiiyixo))^0, 

that, upon observing that 

n-1 

Dp{iy{xo)) = ^ {Dp{iy{xa)), ei) e; + p{iy{xQ)) e„. 



1=1 



becomes 



{Dp{i/{xq)), ei) D'^dn{xa) e; + p{v{xo)) D'^dn{xo) e„ = 0. 



1=1 
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Then 



(4.25) 

for every i — 1, . 



1=1 



+ p{iy{xQ)) {D^dnixa)en, e^) = 0, 



,,n. Hence ( |4.22D follows from ( |4.25| ) and ( |4.2l| ), while ( [4.23D 
follows from ([4.25D and ([4.22D. □ 



Remark 4.19. Formula (4.21) can be written as 

1 



p(t^(a;o)) 



{Dci, ej) , j = 1, . . . ,n - 1, 



where D e Al„_i is the matrix defined in ( 4.15 ). 

5. p-CURVATURES 

Thanks to Theorem 4.16, the function df^ is of class on a neighborhood of 
dQ. We can then define the matrix- valued function 

(5.1) W{xo) = -D^p{Ddnixo)) D^dn{xo) , xo e dfl, 

where D'^dn{xo) = D^dfi{xo). For any xq G dil let T^g denote the tan gent s pace 
to dfl at Xo- Notice that, if w e R", then W{xq)v G T^g, since from ( 4.20 ) and 
Theorem |T](i),(iv) 

{W{xq)v, iy{xo)) = - {D^dn{xo)v, D^p{Ddn{xo j} i^{xq)) 

= -p{iy{xo)) {D'^dn{xo)v, D'^p{v{xo)) v{:xo)) = 0. 

Hence, we can define the map 

(5.2) W{xo):T,g^T,g, W{xo) w ^ W{xo) w, w^T^,, 

that can be identified with a linear application from M"^^ to K"^^. In our set- 
ting, the function W plays the role of the Weingarten map (see Example ^.6| and 
Remark |5.9| below). 

Lemma 5.1. Let x^ G 9f2. Then, in the principal coordinate sys tem at Xq, 
W{xq) — RD where D and R are the symmetric matrices defined in ( 4-1^ ), one 
has 

det[I -tW{xo)] = det[/„_i -tW{xo)], t G M. 
Furthermore, both determinants are strictly positive for t G (— ^^(a;o): ^(^^o))- 

Proof. Consider the principal coordinate system at a;o- Recall that, from the posi- 
tive (— l)-homogeneity of D^p and ( 4.20| ), 

D^p{Ddn(xo))= p{i^)D^p{i^), 

where — ^{xo)- Since D'^p{v) e„ = D"^ p{v) v ~Q (see (p^)), we get 

n-l 

(5.3) {W{xa) e-j, Ci) = -p(t^) ^ {D'^p{v) e;, e^) (D'^dnixo) e^, e;) , 

1=1 



for every i,j — l,...,n — l. From Lemma 4.18 we obtain 

{W{xo)ej, Ci) = {RD),,, 



i,j = l,...,n-l, 
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where D and R are the symmetric matrices in A^„_i defined in (4.15). Since T^^ 
is spanned by (ei, . . . , e„_i), we have that 



W{xo) ^RD. 



(5.4) 

Furthermore 

{W{xo)ek, I') 
hence 

det[I -tW{xQ)] =det[/„_i - tW{xo)] = det(/„_i - tRD), 
and the conclusion foUows from Lemmas 4.1Cl| , 4.11 and Theorem 4.13 . 



□ 



Remark 5.2. Let xq G dfl. From (4.24) we have that 

W{xo)Dp{i^{xo)) = ~D^p{Ddn{xo))D^dn{xo)Dp{i^{xo)) 







that is Dp{v{xt^)) (and hence Dp{DdQ{xo)), by ( 4.20 )) is an eigenvector of W{xo) 
with corresponding eigenvalue zero. Since (Dp{i'{xo)), i'^xq)) = p{i'{xo)) ^ 0, 
then DpjjAxo)) ^ T^o is not an eigenvector for W{xq). On the other hand, from 
Lemma 5.1 we deduce that a number k ^ is an eigenvalue of W{xf)) if and only 
if it is an eigenvalue of W{xq). 

Remark 5.3. Although the matrix W{xq) is not in general symmetric, its eigenvalues 
are real numbers (and so its eigenvectors are real) . This property easily follows from 
(5.4), and the fact that R and D are symmetric matrices. Moreover, since W{x) 
is a continuous matrix- valued function on 9r2, its eigenvalues are continuous real 
functions on dU,. 

Lemma 5.4. Let x G f2 and let xq G n(a::). If k is an eigenvalue ofW(xo), then 
K (x) < 1 . Furthermore, if x €z ^\TZ, then k d^ (x) < 1 . 

Proof. If K < 0, then there is nothing to prove. Assume that k > 0. Then k is an 
eigenvalue of W(xq) if and only if 



det [k In 



W{xo)] = det 



In-l 



-Wixo) 

K 



0. 



From Lemma t3.l| this can happen only if 1/k > l{xo), that is, if kI{xo) < 1. Since 
l{^o) > dQ{x), this proves that Kdn{x) < 1. Furthermore, if a; G il\TZ then 
l{xo) > dQ{x), so that Kdn{x) < 1. □ 

The eigenvalues of W{xo) have a significant geometric interpretation (see Re- 
marks 5.7 and tS.Sl). 



Definition 5.5 (p-curvatures). Let xq G dil. The principal p-curvatures of dfl at 
xq, with respect to the Minkowski norm do, are the eigenvalues ki(xo) < • ■ ■ < 
K„_i(a;o) of W{xo). The corresponding eigenvectors are the principal p-directions 
of dil at Xq. 

Example 5.6 (Euclidean distance). Let K ~ Bi{Q), so that p{S^) ~ |^|. In this 
case, dfi coincides with the Euclidean distance d^ from dil. Let xq G Oft and 
consider the principal coordinate system at xq. Since 



DM) - 1 



1 



1 
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from Lemma 4.18 we recover W{xo) — diag(Ki, . . . , 0) = —D'^d^{xo) and 

W(xa) — diag(Ki, . . . ,k„_i). Hence the principal p-curvatures and the principal 
p-directions correspond respectively to the principal curvatures and the principal 
directions of dfl at xo- 



Remark 5.7 (Normal curvatures). Definition 5.5 is motivated by the following con- 
struction. Let Xo S dQ and consider the principal coordinate system at xo- Let 
X : U M" be a standard parametrization of dfl in a neighborhood of xq satisfying 
(4.11). Let V e M"~^, V 0, and consider the curve on dil 

x{t) = x{tv), t e /, 

where / = {—to, to) and to > is chosen such that tv G Z// for every t Q I. Let 
us denote by V the plane generated by v and v — v{xo) — e„. For r > 0, let 
Zr = Xq + r Dp{v) and let Kr = Zr — r if", so that xq S Kr H dVt. We would like to 
determine r > in such a way that the section H V has a contact of order two 
with the curve x{t) at xq. Such a value of r = r{v) will be the radius of curvature 
of the curve x(t) at xo with respect to the Minkowski distance, hence k{v) — \/r{v) 
will be the curvature of x{t) at xq, that is, the normal curvature of dO, at xq in the 
direction of v with respect to the Minkowski distance. 

For t ^ I, the distance (with respect to the Minkowski norm) from x{t) to the 
section Kr fl V can be estimated by \p'^{zr — x{t)) — r\. The condition of a contact 
of second order becomes 

0= f^pO(zr-x{t)) 

(5.5) rfi' 

= {D^p'\zr ~ xo) i(0), i(0)) + {Dp°izr ~ Xo) , i(0)) . 

Now, Zj. — Xo = rDp(y), so that from homogeneity and Lemma ^.2| we get 
D'p\zr -xo) = - D''p\Bp(v)), Dp\zr - Xo) = Dp\Dp{y)) 



Piy) 



On the other hand, using (L8) and ( 4.11 ) we get 



ji-i 



i(0) = (z;,0), x(0)= 0,^ 



hence (5^) becomes 
(5.6) ^ 



n-1 



{D^p°{Dp(iy))e,, e,) v,Vj 



Pi-) U 



ri-l 



Recalling the definitions of the matrices H and D in (4.16) and (4.15) respectively, 
we conclude that 

- {Hv, v)~^ {Dv, v)^0. 
r p{v) 

Let us define k = k{v) = l/r. The identity above can we rewritten as 
(5.7) R{v) {p{u)Hv, v) - {Dv, v) ^0. 

By construction, the quantity k{v) is the normal curvature of dfi at xo along the 
direction v, with respect to the Minkowski norm. In analogy of what happens in 
the Riemannian case, the principal p-curvatures can be defined as the invariants of 
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the pair of quadratic forms D and p{i') H, that is, the numbers A £ M that satisfy 
the characteristic equation 

det(L> - \p{v)H) =0. 
On the oth er han d, from ( [4.17| ) we have that {p{v) H)^^ = R, where R is the matrix 



defined in (4.15), hence the invariants are the solutions of det(i?L' — A/„_i) = 0. 



Recalhng (5.4), these are the solutions to 

det(W(xo)-A/„_i) = 0, 

that is, the eigenvalues of W(xo). Moreover, if k.i{xo) < ■ • ■ < K„_i(a;o) are the 
eigenvalues of W{xo) and vi, . . . , Vn~i are the corresponding eigenvectors, then it 
is straightforward to check that iii{xQ) = k{vi) for every i — 1, . . . ,n — 1. 

Remark 5.8. Let xq G dQ and let ki^xq) < ••• < Hn-iixo) be the principal p- 
curvatures of dil at .to, that is, the eigenvalues of W{xo)- We claim that 

(5.8) Ki{xo) < k{v) < Kn^iixo) \fv £ Ux„ --^ {w e Tj;gil; ||w|| = l}, 

where ii{v) is the normal curvature of dfl at Xq along the direction v, defined in 
( ^.7| ), and Uxq is the unit tangent bundle of dfl at xq. More precisely, we claim 
that 

(5.9) iiiixo) = min k(v), K„_i(a;o) — max k,{v) . 

Namely, in the principal coordinate system at Xq, k: R"~^\{0} ^ R is a continuous 
function on 5"^^, hence it admits maximum and minimum on S*""^. Let v G S*""^ 



be a maximum point. Since k is positively 0- homogeneous (see (5.7)), then v is also 



a maximum point of k in the open set M" \ {0}, so that Dk{v) = 0. Differentiating 



(5.7) and plugging this identity we obtain 

k{v) p{v) H V — D V — Q , 

that is, ABt{k{v) p{v) H — D) — Q. Thus k{v) is an invariant of the pair of quadratic 
forms D and p{v)H, and hence an eigenvalue of W{xo). Finally, being k{v) the 
maximum of k on 5*""^, we conclude that klv) — k„_i(xo). Reasoning as above, if 
ui is a minimum point of k on S*""^, we deduce that k{w) — ki{xa). 

Remark 5.9 (Anisotropic Weingarten map). We can also give another interpreta- 
tion of p-curvatures (see [|[ §2.3]). We recall that, by definition, the principal 
curvatures of a smooth manifold M at a given point x G M are the eigenvalues 
of the Weingarten map = —dv^ : T^M T^M, where i'{x) is the normal to 
M at X. In out setting, the inward "normal" to Sfi at x is given by the vector 



p{x) :— Dp{v{x)) (see Proposition 4.4 and Remark |4~5| ). Hence the Weingarten 
map, in the Minkowskian setting, should be L^^ = —dp^ '■ Txfl — *■ Txfl. It is readily 
seen that L^^ coincides with the map W{x) defined in (^^). Hence, the principal 
curvatures at x with respect to the Minkowski distance are the eigenvalues of the 
anisotropic Weingarten map W{x). 

Remark 5.10 (Curvatures in Finsler geometry). If is of class C°°, then the 
principal p-curvatures coincide with the notion of curvature in the setting of Finsler 
geometry as defined in 
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6. Regularity of E 

In this section we investigate the structure of the set S. We shall prove that 
S = 7?., and that is a function of class in the whole \ S. Finally we propose 
a regularity result for the set S. The main tool for the proof of these results will 
be a complete description of the set E \ S. 

Definition 6.1 (Focal points). Let us define the map 

(6.1) <S>{x,t) X + tDp{iy{x)), xedn,te'R. 

Fixed xq £ dfl, we say that z = xq + to Dp{iy(xo)) is a focal point of a;o if the map 
$ is singular at {xo,to), that is det £>$(xo, to) = 0. 

In the setting of Riemannian geometry, the point z in the above definition would 
be called a focal point of dft along the geodesic t ^{xo,t), t > 0, normal to dft 
at Xq (see e.g. p^). 

Remark 6.2. Let Y -.U R", U C M"^^ open, be a local parametrization of dVl in 
a neighborhood of xg, such that Y{Q) = xq. The fact that z = xq + t^Dp{v{xQ)) 
is a focal point amounts to say that the differential of the map 5'(?/,t) := Y{y) + 
t Dp{ v{ Y{y)) ) at ( 0, tg ) has not maximal rank, that is det £'\['(0, to) = 0. Recalling 
(|j|), (H.14|) and (jsj), we have 

det £'^'(0, to) = det(/„_i - to W{xq)) = 0, 

hence the focal points of xq G dO. are '^{xq, l/ki{xQ)), if ki{xQ) ^0,i = l,...,n — 1. 
Moreover, by Theorem 4.13| , 1{xq) < l/K„_i(a;o). 

Definition 6.3 (Optimal focal points). Let xq e dft. If k„_i(xo) > 0, we call 
zq = xo + l/Kn-i(a;o) Dp{v{x^)) the first focal point of xo. If, in addition, zq G E, 
we call zo the optimal focal point of xo- We denote by F the set of optimal focal 
points of the boundary points of fJ. 

In the following proposition we prove that E = E U F. 

Proposition 6.4. //a; G E \ E, then kn-i{xQ) dii{x) = 1, where xo is the unique 
projection of x on dfl. 

Proof. Let a; e E \ E, and let (xk) C E be a sequence converging to x. For every 
k £ N, let Uk, Zk £ dfl be two distinct points in 11(2:^), so that 

(6.2) Xk^Vk + dn{xk) Dp{Ddn{vk)) = Zk + dn{xk) Dp{Ddn{zk)) , 
and let 

^ Uk- Zk 



\yk - zk\ 

Since € S*""^ for every fc, we can extract a subsequence (which we do not relabel) 
converging to a point ^ € S^~^. From Lemma ^.41 we have that {yk)k and {zk)k 
converge to xq. From (6.2) we get 



\Vk-zk\ \yk-zk\ 

Since dn is of class in a neighborhood of dVL and p is of class in R" \ 
{0}, passing to the Hmit we get ^ + dQ{x) p{Ddn{xo)) D'^dQ{xo) = 0, that 
is dn{x)W{xo) ^ — ^. Hence l/dn{x) is a nonzero eigenvalue of W{xo). From 
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Remark 5.2, we conclude that l/dn{x) is an eigenvalue of W{xq)^ that is, l/di^{x) = 
kj (xq ) for some je{l,...,n — 1}. 

Finally, let us prove that iij^xo) = kn-iixo)- Assume, by contradiction, that 
Kj{xo) < Kn-i{xo)- Then l/«;„_i(a;o) < l/iij{xo) = dn{xo) < l{xo), a contradic- 



tion (see Remark 3.2). □ 



Proposition 6.5. Let a; G 51 and xq G n(a;). Then E n {xo,x) 



Proof. By Proposition 3.3(iii) we already know that S n (xq.x) — 0. Moreover, 



from Lemma p.4| we have that kn-i{xo) dfi{x) < 1, hence the strict inequality 



Kn-i{xo) dn{z) < 1 holds for every z G {xo,x). This implies, by Proposition 6.4 



that (S \ S) n (a:o, a;) = 0, concluding the proof. □ 
Corollary 6.6. For every xq G dQ, 

(6.3) l{xo) = min{i > 0; xq + tDp{iy{xo)) G S} . 

Proof. Let us denote by X{xo) the right-hand side of (|6.3|). Notice that the point 



zo — Xq + l{xo)Dp{v{xo)) belongs to the ridge set, then, by Remark 4.8, zq G S. 
Hence we have A(a;o) < ^(a^o)- Assume by contradiction that the strict inequality 
holds. Then, if we set z{t) :— xo+t Dp{v{xo)), fori G {\{xo),1(xq)) we would have 
Xq G n(z(i)) and z{\{xq)) G {xq, z(i))nE, in contradiction with Proposition |6. 4 D 



Theorem 6.7. The function I, defined in ( 3.1\ ), and the function t: fl 
defined by 



(6.4) t{x) 




tDp{Ddnix)) eT,}, ifxen\Y,, 
if X eT, , 



(and called normal distance to cut locus of x) are continuous in Q. 

Proof. It is easy to check that, ii x E Q and a;o G n(x), then t{x) — t{xq) — dfi{x). 



Since, by Corollary 6.6, t{xo) — l{xo) = l{x), we get t{x) = l{x) — dnix) for every 



X G f2. From Proposition 3.1C the function I is upper semicontinuous in Q, hence 
the upper semicontinuity of r follows from the continuity of dn. On the other 
hand, the lower semicontinuity of r at a point x G S is trivial. Hence it remains 
to prove that r is lower semicontinuous in x G il \ S. Let {xk)k be a sequence 
in n converging to x. We can assume that Xk G \ S, hence, by definition of r, 
Xk + T{xk) Dp{DdQ{xk)) 6 S. If To := liminffc r(a;i;), passing to the limit we get 
X + tq Dp{Ddfi{x)) G S, that is, r(x) < tq. □ 

Corollary 6.8. The ridge set TZ is closed and 7?. = S. 



Proof. We recall that TZ = {x € il; l{x) = dQ{x)}. Since, by Lemma TZ C ft, 
the fact that TZ is closed follows from the continuity of I and do. The equality 
TZ = Yj now follows from the inclusions S C 7^ C S. □ 

Corollary 6.9. The set S has vanishing Lebesgue measure. 



Proof. It follows from Corollaries 4.15 and |6.8|. □ 



We can now prove a regularity result for d^, which extends Theorem 4.16 to 
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Theorem 6.10 (Regularity of dn in fl \ E). The function dn is of class in 
n\ E. Ifn is of class C^'"" and p € C'='"(M" \ {0}), for some k > 2 and a e [0, 1], 
then dfi is of class C^^"' in $1 \ E. 

Proof. It is enough to prove the resuh in \ S, since the regularity of dn near dVl 
was already proved in Theorem 4.16| . 



Fixed 00 € \ E, we have to show that c?n is of class in a neighborhood of zq- 
Let n(zo) = {xa}, and let F : W — > M", U C R""^ open, be a local p aram etrization 
of dVl in a neighborhood of a;o. Let ^ be the map defined in ( |4.12| ), and let 



to — dn{zo)- By Cor ollary |6.8|, zq does not belong to TZ, that is dn{zo) < l{za) 



l{xo). From Theorem i.l2 we have that det £'5'(xo, to) > 0, hence from the inverse 



mapping theorem it follows that there exists a neighborhood V of Zq such that the 



map y — y{x) is of class C^{V). On the other hand, from Lemma 4.3, Ddn{x) = 
N{y{x))/ p{N{y{x))), x £ V, where N = i' o Y, hence we conclude that Ddu G 
C^V), that is dn e C^{V). 

The last part of the proposition follows from the fact that, if fl is of class C'''" and 
p G C'='"(R" \ {0}), then * G C''-^^°'{U x M). □ 

In literature the set 

Cut{n) = {^{xoJixo)) = xo + lixa) Dp{iy{xa)), xq G dfl} 



is called the cut locus of Q. From (p. 3D, Corollary |6.8|, and Proposition |6^ it follows 
that 

Cut{n) = 7e = E = Eur. 

The regularity of the singular set E and of the cut locus E has been extensively 
studied. For what concerns E, Alberti |l| has proved that it is C^-rectifiable. We 
recall that a subset of R" is C'^-rectifiable, /j G N, if it can be covered by a countable 
family of embedded C*^ manifolds of dimension n — 1, with the exception of a set 
of vanishing 7i"~^ measure. On the other hand, in general the cut locus E is not 
rectifiable. Namely, Mantegazza and Mennucci jl^ have exhibited a set 17 C 
of class C^'^ such that the singular set E corresponding to the Euclidean distance 
d^ has closure E with positive Lebesgue measure. In the same paper the authors 
have proven that, if is an open subset of class C", with r > 3, of a smooth, 
connected and complete Riemannian manifold without boundary, then Cut (SI) is 
(r — 2)-rectifiablc. In the same setting, Itoh and Tanaka have proven that, if 
is of class C°°, then the function / is Lipschitz continuous on d^l. Recently, Li and 
Nirenberg have refined this result, extending it to the case of Finsler manifolds. 
Adapting their result to our setting, they have proved that, if is of class C^'^ and 
p is of class C°°, then / is Lipschitz continuous on 9^7 (see jl^, Theorem 1.5). As a 
straightforward consequence of the fact that / is Lipschitz continuous on dfl, one has 
that the cut locus E has finite 7i"~^ measure and it is C^-rectifiable. Namely, if L 
denotes the Lipschitz constant of I in dfl, we have 7i"~^(E) < L7i"~^(9J7) < +oo. 
Other results in this direction are proved in ||l8|| , again in the setting of Finsler 
manifolds. We mention also |Q, ^ for further regularity results for Hamilton- Jacobi 
equations, and [p^ for rectifiability results of E from an optimal control theory 
viewpoint. 

In the following theorem we shall prove a C°'"-rectifiability result of the cut locus 
in the Minkowskian setting. For what concerns the regularity of the cut locus, our 
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result is finer than the ones given in [Tsj, that, on the other hand, deal with a 
general Finslcr manifold. 

Theorem 6.11 (Regularity of E). Assume that, for some a G [0, 1], 57 is of class 
C^'", and p^p^ G C^'"(R" \ {0}). Then S can be covered by countably many 
graphs of functions of class C'''"(R"~^), with the exception of a set of vanishing 
7i"~^ measure. As a consequence, dimniY,) < n — a, where dimui^) denotes the 
HausdorfJ dimension ofE. 

Proof. Since do is semiconcave in J7, by Theorem 1 in we have that T, is C^- 
rectifiable. Hence it is enough to prove that the set of optimal focal points T can 
be covered by countably many graphs of functions of class C°'", defined on open 
subsets of R"^^ (C°'"-rectifiable for short). Since F C 7?. is contained in the set 

Fq = I xq ^ — - Dp{v{xQ)); Xq e dVl such that K„_i(a;o) > 1/L 

L K„_i(xo) 

of all first focal points with distance from d^l not exceeding the quantity L := 
max{Z(y); y S 957}, then it is enough to prove that Fq is C°'"-rectifiable. Let us 
denote by Sq the set 

So {a;o £ dVL; k„_i(xo) > !/£}. 

Let Ui, . . . , Un C R"^^ be open sets, and Yk'.Uk ^ R", k = 1, . . . ,N, maps of class 
C^'" that parameterize dil. For every k — 1, . . . , N,lct Uk := {y eUk] Yk{y) G Sq}. 
Then^o =Uf=i^fc(f^fc), and 

where, for every k = I, . . . , N, 

■9kiy,t) ■.= Yk{y)+tDp{v{Yk{y))) iy,t)&UkxR. 

Let us fix k ~ 1, . . . , N . We claim that the map K := Hn-i ° Yk is of class C°'" 
in Uk- For every y e Uk let us denote by D{y) and H{y) the matrices defined 
respectively in (4.15) and (4.16) in the point Xq = Yk{y). We recall that D{y) and 
H{y) are the restrictions respectively of D'^d^{Yk{y)) and p° {D p{v{Yk{y)))) to 
the tangent space TY^(y)fl. From Theorem |6.10 , these functions are both of class 
C°'"{Uk). Thus, the matrix 

Hoiy):=p{iy{Yk{y)))H{y), y eUk, 



is of class C'^'"{Uk). Furthermore, from Theorem 2.1 (iv), the minimum eigenvalue 
of H{y) is bounded from below by the minimum for ^ G dK of the radii of curvature 
of dK at ^ (which is strictly positive by the assumption K e C^). From ( ^.l|) we 
conclude that there exists a constant c > such that 

(6.6) {Hoiy)v,v)>c\\v\\^, ^y e Uk, v e W'^'K 

From (pM, K{y) is characterized by 



(6.7) K{y)^ niax J^M!!!^, y ^Uk. 
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Now, let y,y' G^Uk, and denote by vo,Vq G M"^^ two vectors that realize the 
maximum in ( |6.7|) for K{y) and K{y') respectively. Since D,Ho G C'^'°'{Uk), from 
(3.£) we have that 



Kiy) - K{y') < 



{D{y) vo,vo) {D{y')va,vo) 



< C 



{Hoiy) vo, Vq) {Ho{y') vq, Vq) 
ll^o(2/')ll + IIW)ll 



{Ho{y)vo, Vo) {HQ{y')vo, vq) 



y~y'\' 



CS , 

< — \y 



V 



where S ;= sup{||_ffo(^)|| + II^C-^)!!; z G Ut}- Exchanging the role of y and y\ we 
obtain K{y') ~ K{y) < (CS/c^) \y - ?/'|°, hence K G C°-"{Uk). 
Since ^'fc G C^ilAk x M), with bounded derivatives in Uk x [0, L] (see Lemma 4.10 ), 
it is Lipschitz continuous in Uk x [0, L], say of rank L^, hence 



'{(yMK{y))- y^Uk}. 
positive function of class C°'" 

0,a 



7^"-"{vl/fc(y, llK{y))- y G Uk] < W""" 
On the other hand, being K:Uk C M""^ 

on Uk, with K{y) > 1/L for every y G Uk, the function 1/K is of class C"'" on 
Uk- Thus the Hausdorff (n — a)-dimensional measure of the graph of l/ K on Uk is 
finite. Since this fact holds for every k = 1, . . . , iV, from the definition ( |6.5|) of Fq 
we conclude that 7i"~"(ro) is finite. □ 

7. An APPLICATION TO PDEs 

This section is devoted to an application of the previous results to the anal- 
ysis of PDEs arising from optimal transportation theory and shape optimization 
(see |6|| ) ■ The main tool is the following change of variables theorem, based on 



Lemma |4.10| , Theorem [4.13| , and the fact that E has vanishing Lebesgue measure. 
Theorem 7.1 (Change of variables). For every h G L^{il) 

h{x) dx 



(7.1) 



p{iy{x)) 



l{x) 



h{<^{x, t)) det(I„_i - t W(x)) dt 



dK'-^ix) 



/ p{v{x)) / h{'^{x,t))T\{l-tk,{x))dt 
Jon Jo ,_i 



dH^-Hx) 



where $: dfl x 
dn X R. 



is the map defined by ^{x,t) — x + t Dp{v{x)), {x,t) G 



Proof. Let Yk : Uk ^ M", Uk C 



open, k = 1, 



, TV, be local parameterizations 
,iV, let 



of dfl of class C^, such that Ufc=i ^fe(^fc) ~ ^^^r every fc = 1, . 

^-fc : Wfc X M ^ R" be the map 

(7.2) ^kiy,t)^Yk{y)+tDp{iy{Yk{y))), {y,t)eUkxR. 

From Lemma 4.10| and Theorem 4.13 we know that, for every k = 1, . . . , TV, ^E'^ G 
C^Uk X M), and 



det D^fk{y,t) ^ p{:^{Yk{y)) J gf^iy) det[In-i - tWiYk{y))] > 



DISTANCE FUNCTION 



31 



for every y ^Uk and every t £ [0, Z(Yfe(y))), where 
/(y) = det(.gf,(y)), gUy) 



hj = 1, ■ 



, n — 1. 



Let pi, . . . ,pn G dCl ^ M be a partition of unity of dfl subordinate to Yi, . . . , Y^r, 
that is, for every fc = 1, . . . , A^, Pk{x) > for every x € dQ, pk has compact support 



contained in Yk{Uk), Pk °Yk G C'^{iUk), and J2k=iPk{^) = 1 ^'^^ every x G 9ri. 
Moreover, for every /c = 1, . . . , iV let us define the function 

qk{x) 



Pkixo), if xo e Yk{Uk) and x = <i>(xo,t) for some t G [0,;(a;o)), 
0, otherwise. 



and the sets 



Ak^{{y,t)- yeUk, te{0,l{Yk{y)))}, nk = ^k{Ak). 

By construction, it is readily seen that qk = outside Qk for every k, and that 
^k=i qk{x) = 1 for every x G f2 \ S, hence, from Corollary 6.9, for almost every 
a; G fj. The first equality in ( |7.l[ ) is obtained using, for every k = 1, . . . , iV, the 
change of variables '^k on Ak, Fubini's theorem, and the area formula (see [|2[ 
§3.3.3 and §3.3.4]) as follows: 

» N „ ^ r 

/ h{x)dx = ''^^ / qk{x)h{x) dx — ''^^ / qk{x)h(x) dx 

J k—l k—1^^ 

N 

V / qk{'<i^k{y,t))h{-9kiy,t)) detD^kiy,t)dtdy 



fc=i 

N 

E 

fe=i 

E, 



qk{Yk{y))h{^kiy,t))^^^^^Zm^dt 



qk{x) 



l(x) 



h{^{x,t)) p{v{x)) D{x,t) dt 



\Jg^{y)dy 



dW^-^ix) 



h{<^{x,t)) p{v{x)) D{x,t) dt 



dn"^\x) 



where D{x,t) :— det(/„_i — tW{x)). The second equality in (7.1) follows upon 
observing that iii{x), . . . , kn-i{x) are the eigenvalues of W{x) for every x G dfl. □ 



Let us consider the following system of PDEs of Monge-Kantorovich type: 
(7.3) 



' div{v Dp{Du)) = f in f], 
p{Du) < 1 infl, 
p{Du) = 1 in {t; > 0}, 



where the source / > is a continuous function in fJ, complemented with the 
conditions 



(7.4) 



M > 0, > inn, 
u — on dn. 
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The first equation in ( 7.3) has to be understood in the sense of distributions, whereas 
u is a vis cosi ty solution to p{Du) = 1 in the set {v > 0}. We look for a solution 
(w, u) to (7.3)-(7.4) in the class of continuous and non-negative functions. 

Since dfi is a viscosity solution of p{Du) — 1 in Q, u — on dfl, it is clear that 
do is a viscosity solution, vanishing on dfl, of the equation p{Ddii) = 1 in {u > 0}, 
for every fixed continuous function v G C{Q) . Then it is enough to prove that the 
equation 



(7.5) 



div(z;(x) Dp{Ddn{x))) ^ fix) 



in ri, 



has a continuous non-negative solution. More precisely, we are interested in finding 
a function v G (7(17) satisfying 



v{x) {Dp{Ddn{x)), Dip{x)) dx — f{x)(p{x) dx 



(7.6) 



for every ip belonging to the set C^{fl) of functions of class C°°{n) with compact 
support in f2. 



In order to write explicitly a solution v of (7.5), it will be convenient to extend 
the functions <i>, defined on dft in Theorem 



7.1 



and ki, « = 1, 



n — 1, defined on 

dft in Definition 5.5, to f2 \ S in the following way. If x G \ S and n(a;) = {xq}, 
we set 

^{x,t) -.^ xa + (dnix) +t) Dp{iy{xo)), ki{x) -.^ ki{xo) , i = \,...,n-l. 

Since, from Lemma |3.4| the map which associates to every point a; G \ S its 
unique projection on 9il is continuous in \ S, we conclude that the maps Ki, 
i = 1, . . . , n — 1, which are continuous on dVL (see Remark 5.3), are also continuous 
in \ S, and the map $ is continuous in (fi \ E) x K. 

The remaining part of this section will be devoted to the proof of the following 
theorem. 

Theorem 7.2. Let f G C{n), / > 0. Then the function 



(7.7) vf{x) 



/(*(^,t))n 



1 - {dn{x) + t) ki{x) 
1 - di-i{x) ki{x) 



dt ifxeQ\'E, 
if X £ , 



is continuous in Q and sati sfies (7.t) for every if G C^{Q). Here t is the distance 
to cut locus defined in l[6.^). 



Before proving Theorem 7.2, we establish some basic bounds on Vf. 

Lemma 7.3. Let f G C{ri), / > 0. Then, in any set fie := {x G fi; da{x) > e}, 
e > 0, Vf satisfies the bounds 

(7.8) 



0<Vf{x)<\\f\\cin,)Tix)Y[{l + TK^) 

i=l 



Vx G fie , 



where T and are the constants defined by 

(7.9) T = max{r(x); x G dn}, 

(7.10) K ^ — max{[ki{x)] _\ a; G 9^1, « = 1, . . . , n — 1}, 
being [a]_ = max{0, —a} the negative part of a real number a. 
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Proof. By Lemma 5.4 and Corollary 3.8 we have that 1 — {dn{x) +t)Ki{x) > for 
every x G and < i < t(x). Then the function Vf is well defined and Vf > 

in ri. In order to prove (7^) it is enough to observe that 

^-idnix) +t)n {x) ^ 

l-dnix)k,{x) - ^ 
and that x + tDp{Ddi^{x)) G Vl^ for every a; e ile and t £ [0, t{x)]. □ 

Proof of Theorem [7.4 The continuity oi Vf in \ S follows from the continuity 
of the functions k^, i — 1, . . . ,n, in J7 \ S, of the functions / and r in fi, and of 
the function $ i n [f l \ E) x R. On the other hand, the continuity oi Vf on S is a 

consequence of (7^). 

Let if G C^{fl). From the change of variables formula (TJ^) we have that 

f{x)tp(x) dx 



an 



p{v{x)) 



l(x) 



ri-l 



f{^{x,t))^{^{x,t)) Y{{l-tUx))dt 



d7r-^{x). 



Let us compute the term in brackets, integrating by parts, and taking into account 
that (y9($(a;, 0)) — ^{x) — for every x £ dil. Setting $ = <i>(x, t), we get 

n-1 



(7.11) 



/ f{^)^{^)l[{l-tk,{x))dt 

•^0 »=1 

Mx) Mx) 

= / {D^{'^),Dp{v{x))) f{<^{x,s))\\(l-skdx))dsdt. 
Jo Jt 



Noticing that 

Vf{^{x,t)) - 

we obtain 

f{x)ip(x) dx 



x£dn, te [OJix)), 



/ pH / {Dipi^),Dp{iy))vf{^)l[{l-tk,)dt 
Jdn \_Jo 



drt-^{x). 



Finally, again from (7.1), 

f{x)ip{x) dx ^ I Vf{x) {Dp{Ddn{x)), Dip{x)) dx , 
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